Non-mean-field screening by multivalent counterions 
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Screening of a strongly charged macroion by its multivalent counterions can not be described in 
the framework of a mean-field Poisson-Boltzmann (PB) theory because multivalent counterions form 
a strongly correlated liquid (SCL) on the surface of the macroion. It was predicted that a distant 
counterion polarizes the SCL as if it were a metallic surface and creates an electrostatic image. 
The attractive potential energy of the image is the reason why the charge density of counterions 
decreases faster with distance from the charged surface than in PB theory. Using the Monte Carlo 
method to find the equilibrium distribution of counterions around the macroion, we confirm the 
existence of the image potential energy. It is also shown that due to the negative screening length 
of the SCL, —2^, the effective metallic surface is actually above the SCL by |^|. 



I. INTRODUCTION 

In this paper we deal with a problem in which one big 
and strongly charged ion, called a macroion, is screened 
by much smaller but still multivalent counterions, each 
with a large charge Ze (e is the proton charge); for 
brevity, we call them Z-ions. A variety of macroions 
are of importance in chemistry and biology, including 
charged lipid membranes, colloids, DNA, and viruses. 
Multivalent metal ions such as La"*""^, dendrimers, and 
short polyelectrolytes can play the role of the screening 
.Z-ions. 

To illustrate the fundamental aspects of screening we 
use the simple geometry of a solid occupying the half 
space X < 0, whose surface at a; = has a large uniform 
surface charge density —a. The surface charge is screened 
by an aqueous solution of positive, spherical Z-ions with 
radius a, which occupies the rest of space a; > (see 
Fig. [1]). Both the macroion and the aqueous solution 
have dielectric constant e ~ 80. If all of the Z-ions were 
to condense on the macroion's surface, their total charge 
per unit area would equal a. In such a neutral system, 
the concentration of Z-ions N{x) — > at x — > oo. The 
main goal of this paper is to discuss the behavior of N{x). 
The solution of the Poisson-Boltzmann (PB) equation for 
this problem has been known for nearly a century [J, Q ■ 
The Gouy-Chapman solution is 



N{x) 



1 



1 



27rZ2?B {X + x-a)^'- 



(1) 



where A = e/2nalBZ is the Gouy-Chapman length, and 
Ib = e'^/{€kBT) ~ 0.71 nm is the Bjerrum length. We 
have modified the standard Gouy-Chapman formula, tak- 
ing into account the finite radius of the Z-ions, which can 
not approach the surface closer than x = a. 

It was shown d i, i, H, 0, i, i, M, [H that the 
Gouy-Chapman solution fails when both a and Z are 
large enough. The reason it fails is that, in addition to 
A, there is a second length scale in the problem due to 
the discreteness of charge. When the condensed Z-ions 
neutralize the charge of the plane, the two-dimensional 
concentration of Z-ions is n = cr/Ze, and the surface 
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FIG. 1: A stray positive Z-ion (elevated black sphere), at 
a distance x from the surface (thick line) of a negatively 
charged planar macroion (shaded region). Other Z-ions 
(black spheres), condensed at the surface, are on average a 
distance A above the macroion's surface. The dashed line in- 
dicates the average distance, Xc = a + X, from the macroion's 
surface to the adsorbed Z-ions' centers. The stray Z-ion and 
its negative image charge (white sphere) , are equidistant from 
the effective metallic surface, which is shown by the thin line 
at Xfj2et — Xc |^|. 



area per ion, the Wigner-Seitz cell, can be approximated 
as a disc of radius b such that nb'^ = 1/n. Thus, 
b = (7rn)~^/2 = {Ze/nay/^ and 26 is approximately the 
distance between Z-ions. We can construct the dimen- 
sionless ratio 



F = 



ZV/e5 



(2) 



Here F is the dimensionless Coulomb coupling constant, 
or the inverse dimensionless temperature measured in the 
units of a typical interaction energy between Z-ions. For 
example, at Z = 3 and DNA like a — 0.95e nm 
in this paper, we get F = 6.4, A ~ 0.79 nm and b 
nm. Thus, the Coulomb repulsion energy of the Z-ions 
dominates the thermal energy. The result is a strongly 
correlated liquid, which has short-range order similar to a 
Wigner crystal [| i, S i, 0, i, i, M MMMM and is 
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located, practically, at the very surface of the macroion. 
This paper deals only with the strong coupling case: T 3> 
1. Another definition for a Coulomb coupling parameter, 
S — 2r^, was introduced in Ref. and of course, in the 
limit r > 1, S > 1 as well. 

Mean- field treatments, along the lines of PB theory, 
fail at r ^ 1, since when a Z-ion strays away from the 
plane to distances x — a <C b, the electric field of his 
neighbors has no significant x projection. In this range, 
the stray Z-ion is only affected by the electric field of its 
Wigner-Seitz cell (a disc of radius b). Therefore, at x — 
a <^ 6, the surface charge of the macroion is unscreened, 
and the electric field is 2TTa/e. Thus, for < a; — a ^ &, 



N{x) 



ZeX 



exp[— (a; — a)/A]. 



(3) 



(Here following Ref. we used an expression for N{a) 
that ignores the atomic structure of water, while Refs. [3, 
0] tried to take this structure into account). 

Remarkably, the same length A characterizes both this 
exponential decay and the Gouy-Chapman solution, Eq. 
([!]). It is clear that the dramatic difference between the 
exponential decay of Eq. ^ and the power law decay of 
Eq. ^ is due to the effects of correlations. Eq. ^ was 
first obtained in Refs. [1, Then it was re-derived in 
Refs. d, [13] and confirmed by Monte Carlo (MC) sim- 
ulations in Ref. Below we again confirm Eq. ([3]) 
atO<a; — a<C6by MC simulations. However, the fo- 
cus of this paper is on the non-PB behavior of N{x) at 
larger distances x — a > b, which has been predicted in 
Refs. [1, 01 but to our knowledge has never been verified 
analytically or numerically. 

To bring this prediction to mind, let us focus on a sin- 
gle, stray Z-ion located above the macroion's surface at 
X > a + b (Fig. [J). Refs. [1, 0] argue that the nega- 
tive charge of the correlation hole, —Ze, will spread to 
a disc of size x as neighboring Z-ions move to occupy 
the Wigner-Seitz cell the stray Z-ion left behind. This is 
similar to what happens in a metallic surface under the 
influence of an external charge. In fact, this metal-like 
polarization of the SCL on the surface of the macroion 
can actually be described by an image charge that ap- 
pears in the body of the macroion. Because the centers 
of the Z-ions which form the SCL are typically located 
at a distance Xc — a -\- X above the surface (see Fig. [1]) 
it is natural to think that the effective metallic surface 
is at Xmet = Xc and therefore the image is located at 
—X + 2xmet- The attractive interaction energy between 
the stray Z-ion and its image is then [23|, for x—Xmet ?J b, 



Z'e' 



4e(x-Xmet)' 



(4) 



This attractive image interaction, of course, is a correla- 
tion effect. 

The goal of this paper is to verify, by a Monte Carlo 
(MC) simulation and an analytical calculation, that a 
SCL on the insulating surface of a macroion does be- 



have as a metal, and a stray Z-ion has potential energy 
Uim[x). The plan of the paper is as follows. In Sec. |TT] 
we describe our MC procedure. In Sec. IIIH we present 
our MC results for the screening of a spherical macroion 
by Z-ions. To a first approximation they confirm that 
a stray Z-ion at a: > a + 6 has potential energy Uim{x). 
This is the most important result of our paper. 

At a more detailed level, we see in Sec. IIIII that to 
more accurately fit Eq. ([4]) to our MC data the effec- 
tive metallic surface must be lifted slightly above Xc- We 
find that a shift of 0.21 nm provides the best fit. This 
shift is explained in Sec. IIVI where we analytically de- 
rive Eq. ([4]), showing that there is indeed an attractive 
interaction energy between the stray Z-ion and its im- 
age. We further prove that the effective metallic surface 
should be lifted slightly by — ^ = |^|, where 2^ is the lin- 
ear screening length of the SCL. In other words, Xmet, in 
Eq. (HI should be replaced by Xm^t — Xc — 5^ — Xc + 
We show that, theoretically, ^ — —0.20 nm, in reason- 
able agreement with the MC simulation. The fact that 
a Wigner-crystal-like SCL has a negative screening ra- 
dius was predicted theoretically [Tsl and confirmed exper- 
imentally for a low-density two-dimensional electron gas 
in sihcon MOSFETs and GaAs heterojunctions pil. fl7| 
(see also a recent paper [l8|). 

In Sec. |V]we add a small concentration of monovalent 
salt (for example, NaCl) to our system. We show that 
the attractive image interaction persists in this system; 
however, the attraction is weaker due to screening. 



II. MONTE CARLO SIMULATION 

Our setup is similar to the simulations found in 
Refs. [H H [13, HH m, HI. Our system is contained 
within a spherical cell with radius rmax — 10.0 nm. Cen- 
tered within the cell is a spherical macroion with charge 
Qm ~ — 300e and radius Rm = 5.0 nm (—a = —0.95 
e/nm^). The system is populated by 100 Z-ions of charge 
3e and radius a = 0.4 nm. The mobile particles are 
initialized to random non-overlapping coordinates. The 
wall of the spherical cell has a distance of closest ap- 
proach of so that a Z-ion may be placed with its center 
at the wall. Therefore all Z-ions are found at a radial 
distance r within the range Rm -\- a < r < rmax- After 
initializing the system, the total electrostatic energy of 
the system is calculated as 



101 



2e . ^ . dij 



(5) 



where particle i has charge qi {qi = Qm, and for i > 1, 
qi = Ze) located at the center of a hard sphere with 
radius rji (771 — Rm, and for i > 1, rji — a). The distance 
between particles i and j is c?y . The dielectric constant 
is set to e = 80 everywhere and there are no interactions 
with anything outside of the cell. 
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Selecting a particle at random, the MC program at- 
tempts to reposition it randomly within a cubic volume of 
(3.2 nm)^ centered on the particle's current position. The 
total electrostatic energy of the system, 8, is calculated 
after each attempted move. Modeled as hard spheres, 
if any of the particles overlap after an attempted move, 
such that dij < rji + rjj the move is rejected. Additionally, 
any attempted move that places a particle outside of the 
cell, r > Tmax, is also rejected. Otherwise, moves are 
accepted or rejected based on the traditional Metropo- 
lis algorithm. Simulations attempt 52 billion moves, of 
which ~ 4% are accepted, resulting in each particle being 
moved an average of 20 million times. This low accep- 
tance rate is due to most of the Z-ions being condensed 
on the macroion surface where their average separation is 
h = 1.0 nm; one can increase the rate to ^ 8% by shrink- 
ing the volume in which the Z-ion is randomly reposi- 
tioned to (1.6 nm)3. To ensure thermalization, 5 mil- 
lion moves are attempted before beginning the analysis 
of N{r), the Z-ion's radial distribution. 

Following thermalization, N{r) is computed after ev- 
ery 20,000 attempted moves by dividing the simulation 
space around the central macroion into bins that are con- 
centric spherical shells of thickness 0.1 nm, counting the 
ion population within each bin, and then calculating the 
average Z-ion density of each bin. We now introduce the 
empiric mean field potential, (j){r), which corresponds to 
the MC N{r), and is calculated from the radial distribu- 
tions of the ions in the following way. First, the electric 
field is determined at the outer edge of each spherical 
shell by applying Gauss' Law to the integrated charge. 
Then, the potential (j){r) is calculated by discreetly inte- 
grating the electric field in the radial direction. The em- 
piric mean-field potential, 0(r), has nothing to do with 
the PB potential obtained by a solution of the spherical 
PB equation because, due to correlation effects, the MC 
N{r) differs from Eq. ([1]). In the present case, Z-ions are 
strongly condensed at the surface of the macroion, and 
therefore the potential 0(r) decays so fast with r that the 
interaction energy of a stray Z-ion, Ze(t>{r), becomes less 
than ksT already at r > 5.65 nm. 

The main point of this paper is that the concentra- 
tion of Z-ions, N{r), at a distance r from the center of 
the macroion, is only weakly influenced by the empiric 
mean-field potential energy Ze4>{r) and is mostly deter- 
mined by the attractive correlation energy Uc{r). We 
extract Uc{r) — Uc{Tmax) from the simulation data as- 
suming that Z-ions that stray from the macroion surface 
arc Boltzmann distributed according to. 



Ze(f){r) Uc{r) - Uc{rmax) 



N{r) = N{r^ax) exp - 



(6) 

so that the change in the attractive correlation energy 
for a Z-ion moved from r^ax to r is 

Uc{r)-U,{r^ax) = -fcsTln ( J^^] -Ze^{r), (7) 

\N{r.max) J 




Ze 



FIG. 2: The generalization of Fig. [T]to a spherical geometry. 
A stray Z-ion with charge Ze is shown in a cross-sectional 
view at a distance r from the center of a spherical macroion 
with charge Qm, which is covered by condensed Z-ions (black 
spheres). The condensed Z-ions are located at an average 
distance oi Ra = Rm + Xc. = Rm + a + \ from the center of 
the macroion. The stray Z-ion makes a correlation hole with 
charge — Ze, where the concentration of Z-ions is depleted. 
The resulting correlation potential can be modeled as if the 
Z-ion were near a metallic sphere with effective radius Rme.t = 
Rc + Id- The image charges, — Ze and —q located at the 
center, and q' located at a distance r' away from the center, 
are shown by white spheres. 



where we took into account that (t>{rjnax) — because 
our system is neutral. 

We need to recalculate the theoretical form of Uim for 
a spherical macroion geometry (see Fig. [5]), to test that, 
for r - Rmet > b, 



AU{r) = Uc{r)-Uc{rraax) - Uim,ir)-Uimirmax) = 0. 

(8) 

It is known 24] that a charge Ze at a distance r > Rmet 
from the center of a conducting sphere with radius Rmet 
and a net charge of — Ze, induces two image charges 
within the sphere. The charge q' = —ZeRmet/r is lo- 
cated at a distance r' = R'^^^f/r from the sphere's center, 
and the compensating charge —q' is located at the center 
of the sphere. The net charge of the macroion and the 
SCL, — Ze, accounts for the departure of the stray Z-ion 
and is also fixed at the center of the sphere. In the pres- 
ence of these three charges a stray Z-ion, located at r, 
has potential energy given by 



{Zef 



Zeq' 



2{r - r')e 



Zeq' 
2re ■ 



(9) 



The net charge — Ze has fixed magnitude and position 
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because, unlike charges and —q', it is not created by 
the stray Z-ion polarizing the SCL; therefore, the inter- 
action term that involves the net charge does not include 
a factor of 1/2. In the limit x = r — Rm 'ti Rm, we 
recover the planar Ui,n{x) of Eq. (|4]), because Uim{r) is 
dominated by the influence of charge q' 2± —Ze located 

at r' Rrnet - r. 

The first term within the parentheses of Eq. ^ is writ- 
ten for the case when all but one of the mobile charges {Z- 
ions) are located, as in metal, at the surface. This term 
then describes a stray Z-ion's attraction to the fraction 
of Qm left uncompensated due to its departure. In other 
words, this term is used to exclude the stray Z-ion's self 
interaction with its contribution to the mean-field poten- 
tial m. 

To compare Eq. ^ to the simulation in the next sec- 
tion, we take Rmet ~ Rm + Xc = Rc, which aligns the 
metallic surface with the average position of the centers 
of the Z-ions that comprise the SCL (see Fig. [2]). Be- 
cause our macroion is a sphere and not a plane, the 
magnitude of its electric field drops as i? oc l/r^ at 
r > Rm- Therefore, E — 2TTa/e, used to calculate A, 
should be modified slightly since the Z-ion's centers are 
never closer than a to the macroion's surface. We in- 
troduce (Ts = a[RM / [Rm + o)Y to correct the electric 
field. This leads to, = 0.819 e/nm^, = 0.0913 nm, 
Ts — 5.9, and Rmet = 5.49 nm. 



III. RESULTS OF MC SIMULATION 

At/(r) (Eq. ([H])), the difference between the attrac- 
tive correlation energy extracted from the MC simula- 
tion and the result of the image theory is plotted in 
Fig. [3] for Rmet = 5.49 nm (green circles). As expected, 
when r — Rmet ^ i-e. at r < 6.5 nm the difference is 
significantly less than zero since in this range the SCL 
does not function well as a metal due to its discreteness, 
and, therefore, the attractive correlation energy, Uc{r), 
saturates. However, there is also weaker disagreement 
for r > 6.5nm, which decreases with distance from the 
macroion. This suggests that we have improperly iden- 
tified the radius of the effective metallic sphere used to 
calculate Uim{r). To allow for the adjustment of Rmet, 
we introduce the length |^| so that 



Rm.et — Rr 



(10) 



By minimizing the root-mean-square of AU (r) with re- 
spect to 1^1, on the interval 6.4 nm < r < 7.4 nm, we 
determined that |^| ~ 0.21 nm provides the best fit for 
AU{r) = 0. The quality of this fit is illustrated in Fig. 
|3] (red diamonds). This small correction |^| to i?c indi- 
cates that the foundation of Eq. the attractive image 
interaction, is sound. 

In Sec. IIVI we analytically calculate Uim{x) in order 
to find the necessary adjustment in Rmet by considering 
the response of a SCL, made up of adsorbed Z-ions on 
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FIG. 3: The difference, /^U{r) (Eq. between the cor- 

relation attraction energy extracted from the MC simula- 
tion and the result of the image theory, as a function of 
a stray Z-ion's distance from the center of the macroion 
for three different values of the adsorbed Z-ion's screening 
length, 2^. The length, |^|, determines the increased radius, 
Rmet = Rm -I- a -I- a -I- 1^1 , of the effective metallic sphere used 
to calculate Uim{r) (Eq (O). The green circles correspond to 
^ = 0, assumed in the original theory of Refs. 0|, 0|- The red 
diamonds correspond to the best fit to zero, ^ = —0.21 nm. 
The blue squares correspond to, ^ — —0.42 nm and are shown 
for comparison. 
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FIG. 4: Concentration of Z-ions, N{r), as a function of dis- 
tance from the center of the macroion, starting at J?Af + a — 
5.4 nm. The circles represent the data from the MC simu- 
lation. The result of the image theory, Eq. (|lip . is shown 
by short, blue dashes. The medium length, red dashes show 
Eq. p2p . The long, green dashes show the Gouy-Chapman 
solution (Eq. (HJ), with As substituted for A. The error bars 
for the MC data are smaller than the size of the symbols. 



a planar macroion, to the presence of a single stray Z- 
ion above the SCL. It is determined that the SCL screens 
the potential of the stray Z-ion with a negative screening 
length, 2^, where ^ = —0.20 nm. This moves the effective 
metallic surface further away from the macroion's surface 
by 1^1 = 0.20 nm in reasonable agreement with the MC 
data (see Fig. [T]and Fig. [2]). 
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In Fig. m the concentration N{r) obtained from the 
MC simulation is compared to 



N{r) = N{rrnax) exp 



Ze(f){r) Uim{r) - Uira{r max ) 



ksT 

(11) 

which uses ^ = —0.21 nm to calculate Uim(r) (both 0(r) 
and N{rmax) are obtained from the MC simulation). The 
agreement between the MC data and Eq. (jlip is obvious 
when r > 6.5 nm. In Fig. 31 we also compare Eq. ([3]), 
modified for a spherical geometry, 



N{r) 



■ exp 



(r - Rm - a) 



(12) 



to the MC concentration data. At small distances, r — 
Rm + o- ^ b, i.e. r < 5.8 nm, we find good agreement 
with the exponential decay predicted in Refs. [j, 0| and 
confirmed in Refs. d, [l^, Hlf . 

Let us now comment on what happens at larger dis- 
tances from the macroion, which are not shown in Fig. 3] 
and can not be studied well with the small size of the sim- 
ulation cell used in this paper. According to Refs. [3, 
at distances larger than 
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exp 



IMI 



2kBT 



(13) 



from the planar macroion the PB approximation takes 
over and 



N{x) 



1 



1 



27rZ2/B (A 



(14) 



Here fi is the chemical potential of a Z-ion in a SCL. It has 
been shown that for a SCL on a charged background (one- 
component plasma), at 1 < r < 15, ^ is approximated 
well by [ilia], 

fi = -fcBT(1.65r - 2.6iri/4 + 0.26 InF + 1.95), (15) 

where the first term of this expansion corresponds to the 
chemical potential of a Wigner Crystal. For our parame- 
ters, Z = 3 and a = ag = 0.819e nm"^, this leads to the 
length A — 5.18nm. Then, the approximate extension of 
Eq. p4)) to the spherical geometry using X — T — Rm at 
r = 7.55 nm gives \n[N{r)ZeXs/as] = —8.65, very close 
to the MC result -8.77 (see Fig. H]). The idea behind 
the results of Eq. and Eq. is that the correla- 
tion physics at small distances x — a;c ^ A, produces a 
new boundary condition on the concentration of Z-ions 
for the long distance solution of the PB equation 0, Q ■ 

The authors of a recent paper 26] have already studied 
N{r) at large distances by MC simulation in a much a 
larger spherical cell and found that it is in agreement 
with the predictions of the PB approach based on the 
correlation driven boundary condition. They, however, 
did not identify the image domain of distances r which 
we concentrate here upon. Thus, all three asymptotic 



^-2 
Si -4 



^-6 



O MCData 
▲ Eq. (16) 



5.5 



6.0 



6.5 
r (nm) 



°,o,o o,o 
7.0 7.5 



FIG. 5: Concentration of Z-ions, N{r), as a function of dis- 
tance from the center of the macroion. The MC concentration 
(circles; the same data as in Fig. |4]) differs strongly from the 
concentration of Z-ions (red triangles) obtained from Eq. (|16p 
with the empiric mean-field potential potential (j){r). 



regimes, predicted in Refs. [1, namely Eq. ^ a.t x — 
a < 5, Eq. (g]) at & < x - a < A, and Eq. (Hi]) at 
X — a > A are now confirmed by MC simulations. 

We have shown above that the standard mean-field 
theory [U, Q fails to describe screening by multivalent 
counterions. Let us now show that another mean-field 
approximation, which we call the empiric mean-field, fails 
more dramatically. The empiric mean-field potential, 
if){r), was introduced in Sec. |TT] and is obtained using 
the distribution of charge realized in our MC simulation. 
In Fig. [5l we compare the Z-ion concentration obtained 
from the MC simulation to the Z-ion concentration pre- 
dicted using the empiric mean-field potential. 



N{r) = Nq exp 



Ze<f){r) 



(16) 



Here, Nq = 2.18 x 10^^ nm^"^ is the concentration neces- 
sary to normalize the number of Z-ions, in the range 5.4 
nm < r < 10.0 nm, to 100. Clearly, the empiric mean- 
field potential is not self-consistent; Eq. p6|) predicts 
that there are many more Z-ions, at r > 6.0 nm, than 
are actually present in the distribution that produced 
4>{r). The distribution of the Z-ions, for r — Rmet ^ is 
strongly influenced by the attractive correlation interac- 
tion and, therefore, cannot be predicted by the empiric 
mean-field interaction alone. 



IV. THEORY OF IMAGE POTENTIAL AND 
EFFECTIVE METALLIC SURFACE 

In this section we return to the plane geometry of Fig. 
[T] and analytically derive Eq. Q for Uim{x). In the 
process of this derivation, we find the theoretical location, 
Xmet, of the effective metallic surface. The probe charge. 



6 



a stray Z-ion, is positioned far above the plane at x' = 
X ^ b and g = 0, where x' is the axis and g is the radius 
of the cyhndrical coordinate system (a;', g,9). A SCL of 
Z-ions is located in the {g,0) plane at x' = Xc, where the 
typical distance that separates adjacent Z-ions is b. 

The plan of this section consists of, (1) determining 
the analytic solution for the total potential of the system, 
ip{g, x'), (2) presenting it as a sum of two potentials: one 
of the stray Z-ion and the other of the induced charge 
density of the SCL, tlJind{g,x') (the potential of a point 
like image) and (3) finding the position of the effective 
metallic plane, Xmet, so that the attractive interaction 
energy ^ZeipindiO, x) = Uim(x). Below, we show that 
Xmet = Xc — ^, where 2^ is the screening length of the 
SCL, which we also calculate. 

To find the potential ^(g,x') we solve Poisson's equa- 
tion, 



V^^{g,x') 



An 



pig, x'), 



(17) 



where p{g,x') = pext{g,x') + Pind{g,x'), with pext = 
ZeS{g)S{x' — x)/{Trg), and the charge density that is in- 
duced within the SCL is given by 



Pind{g,x') = Ze n^{g,Xc)j-n{0) 5{x' - Xc) 

= Zeip{g,Xc)^S{x' - Xc) 

= -iZef^{g,XcAix' ~xc). (18) 
dp 

Here, n{ip), is the number of Z-ions per unit area as 
a function of the local total potential, ip{g,Xc), and p 
is the chemical potential of the SCL. We consider the 
case, X — Xc b, when the stray Z-ion produces a weak 
potential in the x' — Xc plane {Ze4'{g,Xc)/kBT <C 1). 
This allows us to linearize gind with respect to tl), in Eq. 
([T8| . Rewriting Eq. ((Til) with the help of Eq. pj]) results 
in 



47r 



1 



V'^ii{g,x') = pcxt{g,x') + -il}{g,Xc)5{x' -Xc), (19) 



where, 



dp 



ATr{Ze)^dn 2k 



(20) 



and K is the inverse screening length of the SCL of ad- 
sorbed Z-ions [li,[2|. 

In order to calculate ^ we use p{n) as given by Eq. pS]) 
and the definition of T from Eq. ([2]) and b = (7rn)~^/^. 
For Z = 3, cr = <Ts = 0.819 e/nm^ and e = 80, we find 
that ^ = -0.20 nm. 

In order to solve Eq. (|19p for 4'{g^ x'), we substitute 

/>OC 

i;{g,x')^ / kAk{x')Mkg)dk, (21) 
Jo 



into Eq. (fTg]) . where Ak{x') are the coefficients of the 
expansion and Jo{kg) is the zeroth order Besscl function. 
This yields [28] 



■ip{g,x') 



Ze 



\J {x - x')^ + g^ 



Ze f"" 1 
~h 2fcC + l 



exp \—k{x' -I- a; — 2x^\jQ{kg)dk, 



(22) 



Because the screening length f < 0, the second term 
diverges. To obtain a solution despite this pole, following 
Ref. [l8[, we consider the contribution to ilj{g,x') from 
k <^ only. Such an approach is valid if the stray 

Z-ion, and the observation point large distance 

away from the SCL: [x — Xc), {x' — Xc) ^ This allows 
us to expand l/(2fc^ 1) in Eq. around fc = 0, so 
that 1/(2A:^ -|- 1) ~ 1 — 2k^, and we arrive at 



ip{g,x') 



Ze 



e\/ g'^ + [x - x')^ 
Ze 2{x' + x- 2xc)Zei 



i,J g'^ + [x' + X - 2xc)^ e [^2 -f (x' + X - 2x,)2]3/2 ■ 

(23) 



The first term of Eq. ((23| is the potential created di- 
rectly by the stray Z-ion. The other two terms are the 
first two terms of the expansion of the induced potential, 
4'ind{g,x'), with respect to ^. We are now in a position 
to recast 'i/'md(0, x) at [x — Xc) ^ as being created by 
an image charge a distance s below the stray Z-ion, 



1 



Zeipind{^,x) 



[Zef 



4(a; - Xc)e 
{Zef 
2se 



(Ze)2 



A{x — Xc)'^e 
U^mix), (24) 



where, s — 2{x — Xc + S,)- Specifying that the metallic 
plane must lie half way between the real charge and the 
image charge sets its position at Xmet = x — s/2 = Xc — £,- 
Therefore the effective metallic plane is found a distance 
^ above the plane of the adsorbed Z-ion's centers (Fig. 
[il . This agrees with the statement of Ref. [l^ j that the 
potential created by the stray Z-ion is negative in the 
x' = Xc plane. The theoretical value ^ = —0.20 nm is 
in reasonable agreement with our MC result, ^ = —0.21 
nm (Sec IIIip . Moreover, we have demonstrated that the 
image attraction predicted in Ref. [1, 01 can be derived 
analytically in the limit x ^ b. 



V. SCREENING THE IMAGE BY ADDING 1:1 
SALT 

In this section we modify our system to include a small 
concentration of 1:1 salt molecules such as NaCl. By tak- 
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ing into account the effect of screening on the attractive 
interaction energy between a stray Z-'iow and its image, 
Uirm we obtain a new prediction for the concentration 
of Z-ions, N{r), which is in reasonable agreement with 
the new MC resuhs. In calculating Uim we assume that 
the concentration c of 1:1 salt is so small that the total 
potential ip at any point in the bulk solution (r > 6.0 
nm) obeys the linearized Poisson-Boltzmann equation, 



(25) 



where l/nb is the Debye-Hiickel (DH) screening length, 



1 ^ eksT 

Kb 



(26) 



The exact solution of Eq. ((25)) for a point charge a dis- 
tance r—Rmet away from the surface of a grounded metal- 
lic sphere in a weak electrolyte is known [2^; however, 
we will avoid the complexity of this solution and approx- 
imate the spherical macroion and its SCL as a grounded 
metallic plane. As seen in Fig. [6l where Eq. ([4]) is used to 
calculate N{r) (short blue dashes) for c = 0, one obtains 
reasonable agreement with the MC N{r) (circles) without 
using Eq. ([9]) as we did in Sec. IIIII This demonstrates 
that the influence of the total central charge, —Ze — q', is 
very small. The reason for this is that when a stray Z-ion 
is close to the macroion surface, the total central charge 
is much smaller than the image charge q' . Additionally, 
the central charge is much farther from the stray Z-ion 
than the image charge q' . When the system includes 1:1 
salt, the influence of the total central charge is further 
reduced due to screening. 
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MC Data 
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Eq. (28) 
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FIG. 6: Concentration of .Z-ions, N{r), in a weak electrolyte 
solution, as a function of distance from the center of the 
macroion. The shapes represent the data from the MC sim- 
ulations for two different concentrations of 1:1 salt: 15.2 mM 
(triangles), and (circles). The result of the screened im- 
age theory, Eq. (f28)l with x — r — Rmet + Xmet, is shown 
by medium length red dashes for c=15.2 mM, and short blue 
dashes for c=0. The error bars for the MC data are smaller 
than the size of the symbols. 



Consider a Z-ion which is submerged in a weak elec- 
trolyte solution with dielectric constant e, and is a dis- 
tance X — Xmet above a grounded metallic plane located 
at Xmet (see Fig. [T]). For this system, the solution to Eq. 
P5p . subject to the boundary conditionj^(a;met) = 0, can 
be found using the method of images |2!| . To satisfy the 
boundary condition, the image potential must exactly 
cancel the potential of our Z-ion in the metallic plane. 
Such an image potential is provided by the DH screened 
potential of a charge {—Ze) located at x' = —x + 2xmet- 

The interaction energy of a stray Z-ion with its 
screened image is now readily calculated and is given by. 



Uimi-^) — 



Z^e^ 



4e(x - Xmet) 



exp{-2Kb[x - Xmet])- (27) 



In the limit of infinite dilution c — > , or equivalently 
Kb — > 0, Eq. (|27| is equal to Eq. (|4]), as expected. On- 
sager and Samaras ^30j obtained the same result, but 
with the opposite sign, for an ion's interaction with its 
image at an electrolyte-air interface resulting in a repul- 
sive force. To compare Eq. ((27|) to the spherical geometry 
of our MC simulation, we take x — t ^ Xmet — Rmet using 

Rmet = Rc+ 

The MC simulation described in Sec. |TT] was modified 
to include M 1:1 salt molecules, resulting in M ions of 
charge e and M ions of charge — e. All of the monova- 
lent salt ions have their charge at the center of a hard 
sphere with radius rj = 0.2 nm. We studied a 1:1 salt 
concentrations of 15.2 mM corresponding to the addition 
of M — 34 salt molecules to the solution. The following 
changes were made to the MC simulation to properly in- 
corporate the new ions. The sum used to calculate the 
total electrostatic energy of the system (Eq. ([5])) was 
changed to include the monovalent ions, and the mono- 
valent ions were also incorporated into the calculation of 
the empiric mean- field potential ^(r). 

In Fig. [51 the concentration N(r) obtained from the 
MC simulations is compared to 



N{r) = 7V(ro)exp 



Ze[0(r) - 0(ro)] 

Uimjr) - Uimjro) 
knT 



(28) 



which uses ^ = —0.21 nm to calculate Uim{x) (Eq. (P?]) ) 
(both (j){r) and N{ro) are obtained from the MC simula- 
tions). Because there are no screening particles outside 
of the simulation cell, Z-ions near the wall of the cell 
are repelled from this interface [3l| even though there is 
not a jump in the dielectric constant. To keep this effect 
separate from the image interaction of interest, we chose 
our reference point at rg = 8.05 nm instead of Vmax- Be- 
cause the stray Z-ion's attraction to the image is reduced 
by screening, we see that in Fig. [6] the concentration of 
Z-ions is higher at distances r > 60 nm when 1:1 salt 
is present in the solution. Even with the addition of 1:1 
salt to the solution the agreement between the MC data 
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and Eq. ((28)) for r > 6.5 nm is reasonable, demonstrating 
that the metaUic behavior of the SCL on the macroion 
surface survives, and that the image attraction is still 
important to determining the Z-ion's concentration. 

VI. CONCLUSION 

To summarize, we have studied the role of correla- 
tions among adsorbed Z-ions in attracting stray Z-ions 
and influencing their distribution in the screening atmo- 
sphere. Adsorbed Z-ions on the surface of the macroion 
form a strongly correlated Hquid (SCL). The SCL acts 
as an effective metallic surface for Z-ions that stray from 
the macroion surface to distances larger than the aver- 
age distance between Z-ions of the SCL. Using Monte 
Carlo (MCj_ simulations, we verified the theoretical pre- 
diction [3', '3] that a stray Z-ion is attracted to its elec- 
trostatic image created behind the effective metallic sur- 
face. As a small correction to Refs. [1, 0]) however, our 
MC simulation showed that the effective metallic surface 
is not aligned with the average position of the adsorbed 
Z-ion's centers, but is slightly above the adsorbed Z- 
ion's centers. This offset was calculated analytically to 
be 1^1, where 2^ is the screening length of the SCL. Our 
analytic theory is in reasonable agreement with the MC 
data. Extending the original image theory of Refs. 
we demonstrated that the attractive image interaction, 



while screened, persists in a weak electrolyte. 

In Ref. the attractive image interaction, which we 
have studied here, was used to interpret the origin of 
the negative chemical potential of the condensed Z-ions 
(Eq. [15]). As a stray Z-ion, attracted to its image, ap- 
proaches the surface of the macroion it reaches a dis- 
tance, r ~ Rm + b, where the SCL fails to act as a good 
metal and the correlation attraction Uc{r) saturates at 
/i ~ —Z^e^/Aeb; this saturation can be seen in Fig. |4]as 
the growth of AU{r) at r < 6.5 nm. It is this negative 
chemical potential, brought about by the attractive im- 
age interaction, at r > Rmet + b, which drives charge in- 
version (over compensation of the macroion's bare charge 
with condensed Z-ions), a phenomenon that has gener- 
ated much interest i, i, i, i, 0, H i, [M [HI- Thus, 
we believe that this paper helps to clarify the origin of 
charge inversion. 
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